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Abstract 

The Improved Quantum Molecular Dynamics model is further developed by introducing new 
parameters in interaction potential energy functional based on Skyrme interaction of SkM* and 
SLy series. The properties of ground states of selected nuclei can be reproduced very well. The 
Coulomb barriers for a series of reaction systems are studied and compared with the results of the 
proximity potential. The fusion excitation functions for a series of fusion reactions are calculated 
and the results are in good agreement with experimental data. 
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I. INTRODUCTION 

Recently, the study of the mechanism for heavy-ion fusion reactions at energies near the 
Coulomb barrier, especially, the mechanism for the enhancement of fusion cross sections 
for neutron-rich systems has attracted a lot of attention. The knowledge of the mechanism 
of the enhancement of fusion cross sections for heavy and neutron-rich systems is useful 
in the synthesis of superheavy elements. In heavy-ion fusion reactions, the excitation and 
deformation of projectile and target, neck formation and nucleon transfer strongly influence 
the dynamics of fusion processes. These effects are more pronounced at near and below the 
fusion barrier. For heavy systems, in order to consider these effects a very large number of 
degrees of freedom of motion is involved and the situation becomes very complicated. Thus 
one will meet great difficulty by macroscopic dynamics models U] in which only few degrees 
of freedom of motion are included. The difficulty is also encountered by the fusion coupled 
channel model P] for it is quite difficult even impossible to include such a large number of 
possible channels in practical calculations. Therefore, it is highly requisite to develop a 
microscopic dynamical model suitable for studying heavy ion fusion reactions by which one 
can consistently take account of the dynamical deformation, particle transfer, isospin and 
mass asymmetry effects, etc. In our previous work^, |^, an improved quantum molecular 
dynamics (ImQMD) model was proposed. Main improvements in the ImQMD model are 
as follows: The surface and surface symmetry energy terms are introduced in the potential 
energy part; A system size dependent wave packet width is introduced in order to consider 
the evolution of the wave packet width; An approximate treatment of anti-symmetrization, 
namely, a phase space occupation constraint is adopted . With this model we have studied 
the dynamical evolution of the fusion barrier as well as the development of the neck in fusion 
reactions of ^°'^^Ca+^°'^^Zr. However, more tests of our model are needed. Before that one 
of the most urgent problem has to be solved is to increase the time of keeping an individual 
nucleus to be stable and in a good shape (close to ground state shape) without emission of 
nucleons in order to study the dynamical process for fusion reactions of heavy nuclei. As is 
well know that the time scale for the formation process of a compound system in a fusion 
reaction of heavy nuclei is about thousands fm/c, and also the time scale of quasi-fission 



which reduces the fusion probabihty substantially in fusion processes of heavy nuclei is about 
several thousands fm/c or longer. To meet this requirement, in this paper we develop an 



updated version of improvec 



quantum molecular dynamics model named ImQMD-II based 



on our previously workj^ ^\. In this paper we mainly devote to modify the parameters 
of the potential energy functional in the model. As is well known that the parameters of 
the potential energy functional in the QMD model are obtained based on Skyrme forces. 
There are quite a lot of new versions of Skyrme forces having been developed following the 
development of the knowledge of nuclear equation of state. The newly developed Skyrme 
forces such as the SLy series are designed to study the properties of nuclei away from the 
/3-stability line in addition to nuclei along the /^-stability line. It has been shown that these 
modern Skyrme forces can describe the properties of the nuclei away from /5-stability line 
better than the old Skyrme forces^]. Therefore it seems to us to be worthwhile to try new 
parameters for the ImQMD model based on the modern parametrizations of Skyrme forces 
such as the SLy series as well as the other popular versions. 

The paper is organized as follows. In Sec. II, we briefly introduce the formalism of the 
ImQMD-II model and the improvements. The applications of this model to the fusion 
reactions at energies near the Coulomb barrier are reported in Sec. III. Finally, the summary 
and discussion are given in Sec. IV. 

II. FURTHER DEVELOPMENT OF THE IMPROVED QMD MODEL 

In this section we introduce the updated ImQMD model(ImQMD-II) in more details. 
First, a brief introduction to the ImQMD model is presented. Then, we will give the new 
development of ImQMD. Finally, the calculation results about the properties of selected 
nuclei with the ImQMD-II model are given. 

A. Brief Introduction of the ImQMD model 

For readers convenience, let us first briefly introduce the ImQMD model. In the ImQMD 
model,the same as in the original QMD modelj^, 0, 0, 13, each nucleon is represented by a 
coherent state of a Gaussian wave packet 

'^^W = (^^-P[-^^^ + ■ P.], (1) 
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where, rj, pj, are the centers of i-th wave packet in the coordinate and momentum space, 
respectively, represents the spatial spread of the wave packet. The total N-body wave 
function is assumed to be the direct product of these coherent states. Through a Wigner 
transformation, the one-body phase space distribution function for N-distinguishable parti- 
cles is given by: 

/(r.P)-i:^expi-ii^-^(p-p.n. (2) 

For identical fermions, the effects of the Pauli principle were discussed in a broader 

n 

context by Feldmeier and Schnack llj. The approximate treatment of anti-symmetrization 
is adopted in the ImQMD model by means of the phase space occupation constraint method 
. The density and momentum distribution functions of a system read 

p(r)= / f{r,p)d'p = J2p,{r), (3) 
■' i 

9ip) = J f{r,p)d\ = Y.9iiP), (4) 

i 

respectively, where the sum runs over all particles in the system. Pi{r) and gi{p) are the 
density and momentum distributions of nucleon i: 

A(r)=(^J^exp[-i^:^], (5) 



^^(P) = (^;^^-P[ (6) 

where and ap are the widths of wave packets in coordinate and momentum space, respec- 
tively, and they satisfy the minimum uncertainty relation: 

ar(Tp = ^. (7) 

The propagation of nucleons under the self-consistently generated mean field is governed 
by Hamiltonian equations of motion: 

OH . dH 

The Hamiltonian H consists of the kinetic energy and the effective interaction potential 
energy: 

H = T + U, (9) 



T = y^. (10) 

The effective interaction potential energy includes the nuclear local interaction potential 
energy and the Coulomb interaction potential energy: 

U = Uioc + Ucoul. (11) 

And 

Uioc = J Vioc{r)dr, (12) 

where Vioc{r) is potential energy density. 

The potential energy density V;oc(r) in the ImQMD model reads 



Vioc=^^ + ^C^ + f^(Vp)^ + 9.P^ + ^{p'- ^s{ypf)5\ (13) 
2 po 7 + 1 Po 2po Po 2po 



where 5 = , The first three terms in above expression can be obtained from the poten- 

Pn+Pp ^ ^ 

tial energy functional of Skyrme forces directly. The fifth term is the symmetry potential 
energy part where both the bulk and the surface symmetry energy are included. In addition, 
we introduce an extra small correction term Vr = Qt—t- (named tau term) in the poten- 
tial energy functional. Inserting expression (13) into (12), we obtain the local interaction 
potential energy omitting self-energies: 

Ui^^ = ^yyPii + ^—y(yPii] (14) 



^ i j^i Po i Po / ^ i j^i Po 



where 



and ti = 1 for protons and —1 for neutrons. One should notice that the third term in (14) 
comes from both surface term and the correction to the second term of (13)(seey]), and 
thus go is actually treated as a parameter in the model. 

The Coulomb energy can be written as a sum of the direct and the exchan ge c ontribution, 
and the latter being taken into account in the Slater approximation [^,13, 1^ : 

Ucoui = \j pp{v)y-^pp{v')dvdv' - [lY' I p'J'dR. (17) 



where pp is the density distribution of protons of the system. The colhsion term and phase 
space occupation constraint can also readjust the momenta, but the former plays a very 
small role in low energy heavy ion collisions _jind the latter only happens occasionally. The 
phase space occupation constraint method 5] and the system-size-dependent wave-packet 
width are adopted as that did in the previous version of ImQMD0,Q- 

B. The new development in ImQMD(ImQMD-II) model 

The new development of the ImQMD model mainly are reconsidering the parameters 
in the potential energy functional (see expression (13)). In the expression (13), The first 
three terms can be obtained from the potential energy functional of a standard Skyrme 
interaction directly. The parameters a, /3, 7 and Qsur can be related to the parameters of 
Skyrme interactions by 

= 4^3, (19) 
7 + 1 Po 16 

and 

T^ = i(9^i- -4^2x2). (20) 

Where the parameter 7 is taken the same value as in the Skyrme interaction. The linear 
density dependence of the symmetry energy term is taken and the parameter is fixed by 
the symmetry energy coefficient. In Table. 1 we list the a, /?, and 7 parameters used in 
the QMD modellSlI and the corresponding values obtained from various versions of Skyrme 



interaction 



0,Q,Tii 



17 



18|. From the table one can find the parameters of QMD (hard) is 



close to that of SIII, and the soft one is close to SkP. 

Table 1 



Parameter 


a{MeV) f3{MeV) 


7 


gsuriMeVfm^) 


QMD(hard) 


-124.0 


71.0 


2 




QMD(soft) 


-356.0 


303.0 


7/6 




SIII[15] 


-139.6 


71.4 


2 


20.14 


SkP [18] 


-362.8 


309.6 


7/6 


19.84 


SkM[16j 


-327.3 


258.0 


7/6 


20.32 


SkM*[17] 


-327.3 


258.0 


7/6 


21.82 


SLylOfe] 


-310.2 


228.8 


7/6 


21.57 
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The SIII parametrization was proposed in 1975 and with it one could describe the ground 
state properties of spherical nuclei very well. However, the incompressibility modulus of 
symmetric nuclear matter obtained from SIII (Koo ~ 365MeV) is too high^^. Taking 
this into account, Krivine et al. derived the SkM interaction at 1980. Later on, detailed 
.ud,e. of tKe fiss.on ba..*.Q i„ the ae.n,de .eg., resulted in a .no.e .e«.ed value of 



the nuclear surface tension. Then SkM* [17| was derived which gave a more refined surface 
tension. Nowadays, the symmetry energy part of the interaction attracts a lot interests as 
nuclei further away from the stability line can be expected to be produced with the coming 
radioactive beam facilities. Putting more emphases on the isospin degree of freedom, a 
series of sets of SLy parametrizations were proposed at the late of 1990 's for reproducing 
the properties of the nuclei from the (3 stability line to the drip lines 0]. 

Considering the successes of SkM* in describing the surface tension and SLy in describ- 
ing the properties of nuclear systems far away from f3 stability line, the SkM* and SLy 
parametrization can provide us with a reference to adjust the new ImQMD parameters. 
Concerning the symmetry energy term, the linear density dependence of the bulk symme- 
try energy term is taken and the parameter is fixed by the symmetry energy coefficient. 
The surface symmetry energy term is also introduced. This term is important for having 



a correct neutron skin, which was introduced in the liquid-drop model 



interactions 16 



21 



221 . For Skyrme 



23l |. the surface-symmetry term can be extracted and reads 

Usurf-symm = f [V p{r)]^ 6 {rf dv . (21) 

zpo J 

It modifies the symmetry potential at surface region, and therefore it is especially important 
for correctly describing the neck dynamics in fusion reactions of neutron-rich nuclei. In our 
previous work this term was introduced and it was found that this term played a role in 
the fusion dynamics. We will study this effect in the following section. Taking SkM* and 
SLy parametrizations as the reference we propose a new set of parameters for the ImQMD 
model by reproducing the properties of ground state of selected nuclei ^''^Pb,^°Zr,'^°Ca and 
and the fusion cross sections of ^°Ca+^sCa|2^, ^oCa+^^'^^ZrQ- The new set of ImQMD 
parameters named IQl is listed in Table. 2 

Table 2 



Para. 


a{MeV) (3{MeV) 7 


goiMeVfm^) g^MeV) 


7] 


Cs{MeV) 


K-sifm^) Poifm ^) 


IQl 


-310.0 258.0 7/6 


19.8 9.5 


2/3 


32.0 


0.08 0.165 
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From Table 1 and Table 2, one can see that the parameter set of IQl is generally close 
to SLylO and SkM*. 



Table 3 



Nuclei 


208pb 


i^oCe 


i32Sn 


i^^Sn 


902r 


56Ni 


^^Ca 


^oCa 


160 


Binding 


exp.\17]{MeV) 


7.87 


8.38 


8.35 


8.53 


8.71 


8.64 


8.67 


8.55 


7.98 


Energy 


lQl{MeV) 


7.77 


8.35 


8.27 


8.51 


8.71 


8.63 


8.67 


8.65 


8.23 


RMS 


exp.[17](/m) 


5.50 


4.88 






4.27 


3.75 


3.48 


3.49 


2.73 


Radius 


IQl(/m) 


5.51 


4.87 


4.79 


4.55 


4.25 


3.71 


3.54 


3.44 


2.72 



Fig. 1 

In Table 3 we list the binding energies and the root-mean-square charge radii of ^°^Pb, 
^^°Ce, i32Sn, ii^Sn, ^^Zr, ^^Ni, ^^Ca, ^°Ca, and ^^O calculated by the ImQMD-II model with 
parameter set IQl. The experimental data are also listed for comparison and it is shown that 
the calculated results are in good agreement with experimental data. In Fig.l, we present 
the time evolution of binding energies and root-mean-square charge radii for ^^Zr and ^osp]-, 
calculated by the ImQMD-II model with IQl parameters. One can see that their binding 
energies, and root mean square charge radii remain constants with a very small fluctuation 
and the bound nuclei evolve stably without spurious emission for a period of time of about 
3000fm/ c, which is essential for applications to fusion reactions of heavy nuclei as is discussed 
in the introduction. 

III. Applications to fusion Reactions near the Coulomb Barrier 

In this section, we show the calculation results of Coulomb barriers and fusion excitation 
functions for a series of fusion systems by means of the ImQMD-II model with parameters 
of IQl . 

A. The Coulomb Barrier 

The interaction potential V{R) is defined by 

V{R) = E^2{R) - ^1 - E2. (22) 



Here R is the distance between the centers of mass of projectile and target. Eui^R) is 
the total energy of whole system, while Ei and E2 are the energies of projectile(like) and 
target (like) part, respectively. For kinetic energies, the Thomas- Fermi approximation is 
adopted as mentioned in ref.|2a]. By using the ImQMD model, both the static and dynamic 
Coulomb barrier can be calculated. For the static Coulomb barrier case, the static density 
distribution which is the same as the initial density distribution of projectile and target is 
adopted, while for the dynamic Coulomb barrier case the density distribution of the system 
changes dynamically due to the interaction between the reaction partners. 

We show the static Coulomb barriers calculated by the ImQMD-II model with param- 
eter set of IQl for ^^Ca+'^^Ca, ^^Ca+^^Zr, ^^O+^^O and ^^0+^°^Fh in Fig.2. The results 
calculated from proximity potential 0] are also shown in the figures. One can see that 
the Coulomb barriers calculated with ImQMD-II are in good agreement with those from 
proximity potential js^]. 

Fig. 2 



Fig. 3 

In addition, the effects of the mass asymmetry of projectile and target on the static Coulomb 
barrier are studied through calculating the static Coulomb barriers for ^^^I+^'^^I, ^^Cr+^°^Pb, 
^^S+^^°Th and ^^C+^^°Fm fusion systems which can form the same compound nucleus ^^^Sg. 
In Fig. 3, the solid curves denote the static Coulomb barriers calculated by ImQMD-II with 
parameters IQl and the dashed curves denote the results from the proximity potential. From 
Fig. 3 one can get two points: the first one is that the results from ImQMD-II are in good 
agreement with those from proximity potential when two nuclei do not overlap too much in 
space. The proximity potential may not be able to give an accurate result at the overlapping 
region where the ImQMD model is applicable. Clearly, the results for the overlapping region 
are more interesting, especially for the cases of heavy systems. The second one is that with 
the increase of the mass asymmetry of projectile and target (from i3ij_|_i3ij ^^C+^^^Fm) 
the height of the Coulomb barrier decreases gradually and the capture probability should 
be enhanced consequently. 26^ . 

B. Fusion cross sections 
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In the ImQMD model, we first create certain reaction events for eacli incident energy E 
and impact parameter b (in tis work tlie number is 100) and then counting the number of 
fusion events, we obtain the probabihty of fusion reactions, gfus{E, b), then the cross section 
is calculated by using the expression^]: 

^max 

afusiE)=27r J bgfusiE,b)db = 2nY,bgf^,{E,b)Ab. (23) 



At a certain incident energy, the probability of fusion reactions, gfus{E,b), decreases with 
the increase of impact parameter as shown in Fig.4. This is because the interaction between 
two nuclei decreases gradually from central collisions to peripheral collisions, and following 
it the probability for fusion reactions decreases and that for elastic scattering processes 
increases. In addition, the gfus{E, b) decreases quickly with the decrease of incident energies 
within the energy range interested in this work. Fig. 5 presents the evolution of the fusion 
probability with impact parameter at E^m. = 60, 58, 54, 52Mel^, respectively. One can see 
from the figure that the fusion probability falls when the incident energy decreases from 
60MeV to 52MeV and for energies below the Coulomb barrier, fusion events generally only 
occur at central collisions. 



Fig 


-.4 




Fig 


-.5 




Fig 


-.6 



In Fig. 6, we show the fusion excitation functions for ^°Ca+^^Ti 28|, 46'jj_j_46'jj 



iQ, 



^0Ca+90'962r|23, 32.34g+89Y|2i|, and '^^Si,^^C\+^^ZT at energies near the Coulomb bar- 
rier, and the experimental data are also presented for comparison. In the figure the triangles 
denote the results of ImQMD-II with parameters of IQl and the circles denote the experi- 
mental data and crosses denote the results of one-dimension WKB approximation js^. One 
can see from the figure that calculation results of ImQMD-II for the fusion excitation func- 
tions are generally in good agreement with experimental data, which implies that our model 
is quite reasonable. Now let us study the influence of the surface-symmetry energy term on 
the fusion cross sections, we show the calculation results of the fusion excitation functions of 
^°Ca+^^Ca [2^ at energies near the Coulomb barrier without and with the surface-symmetry 
energy term taken into account in Fig. 7. One can find from the figure that at incident en- 
ergies above the Coulomb barrier the fusion cross sections calculated under two cases are 
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approximately equal, while at incident energies below the barrier the difference between two 
cases become obvious, and for this case the fusion cross sections without surface-symmetry 
energy term taken into account are obviously larger than the experimental data and those 
with surface-symmetry energy term taken into account can reproduce the experimental data 
well. We have found that the N/Z ratio at neck region is enhanced at the early stage for 
neutron-rich fusion processr4], which is driven by the symmetry potential. It lowers the 
fusion barrier and consequently enhances the fusion cross section. If only the bulk term 
is taken in account, the effect of the symmetry energy term becomes too strong at surface 
region and therefore a surface-symmetry energy term should be taken into account in order 
to reduce the effect of the bulk symmetry energy term at surface region. It is especially 
important for having a correct neck dynamics in neutron-rich nuclear fusion reactions. 



FiK.7 



Fig.8 



For further testing the reliability of ImQMD, we calculate the excitation function for 
fusion reactions of the neutron- rich radioactive beam ^^^Sn bombarding on the neutron-rich 
target ®^Ni at energies near the Coulomb barrier, which was recently measured 3^. Fig.8 
shows the comparison of our calculation results and the experimental data as well as the 
results from coupled-channel calculations j;33]. The solid circles denotes the experimental 



data, the triangles denote the results of ImQMD-II with parameters of IQl. The dashed 
and solid curves denote the results of coupled-channel calculations with including inelastic 
excitation (IE) of the projectile and target and IE plus neutron transfer (n&IE), respectively 
(see ref.jSSi] and references therein). From Fig.8 one can find that the coupled- channel 
calculations significantly under-predicted the sub-barrier fusion cross sections 3^ . This is 
because with the coupled-channel model it is difficult to consider all degrees of freedom of 
motion, the number of which is extremely large in sub-barrier fusion reactions for heavier 
systems, while with the ImQMD model, all degrees of freedom of motion are self-consistently 
included. Thus, this model may possibly provide us with a useful approach to explore 
the mechanism of the capture process in the synthesis of superheavy elements. The work 
concerning this aspect is in progress and the results will be reported in the future work. 



Fig.9 
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As a test we also make applications of our model to describe the charge distributions 
of fragments in multifragmentation processes. In Fig. 9 we show the charge distribution of 
fragments by using the ImQMD-II model with parameter set of IQl for Ca+Ca a.t E = 
35MeV/nucleon[34] and Xe+Sn at E = 50MeV/nucleon[35]. One can find that the charge 
distribution of fragments, especially the number of intermediate mass fragments is in good 
agreement with experimental data. It is well known that the number of intermediate mass 
fragments is usually being under-predicted in QMD model calculations. Our results are 
encouraging and it seems to us that the model can be used for heavy ion collisions at both 
low energies and intermediate energies. 

IV. Conclusions and Discussion 

In this work, we have made further improvements in the ImQMD model and proposed 
a new version, namely, the ImQMD-II model in which a new small correction term and a 
surface symmetry energy term are introduced in the potential energy functional in addition 
to the terms adopted in the normal isospin dependent QMD model. A parameter set for 
the potential energy functional based on Skyrme interaction of SkM * and SLy series is 
introduced. By using the new version of the ImQMD model, the ground state properties of 
a series selected nuclei can be described very well. The time evolution of individual nuclei 
can keep stable for about several thousands fm/c which roughly fits the requirement for 
study of fusion reactions of heavy nuclei. We have shown that with this model both the 
Coulomb barriers and the fusion excitation functions for a series of fusion systems (including 
neutron-rich radioactive beam ^^^Sn +®^Ni fusion reactions) at energies near the barrier can 
be reproduced very well. Our study shows that the microscopic dynamical model such as 
the ImQMD model has an advantage of taking account of the dynamical effects such as the 
excitations of projectile and target (deformation and vibration), neck formation, isospin and 
mass asymmetry, etc. simultaneously and thus offers a useful way to study fusion reactions 
of heavy nuclei. Furthermore, we have also shown that this model seems to work well on 
the charge distribution of fragments in multifragmentation processes. 
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CAPTIONS 



Fig.l The time evolution of binding energies and root-mean-square charge radii for ^°Zr 
and ^''^Pb calculated by the ImQMD-II model with parameter set of IQl. 

Fig. 2 The static Coulomb barriers of ^°Ca+^8Ca, ^^Ca+^oZr, ^^O+^^O and ^^^O+^ospb. 
The solid and dashed curves denote the results of ImQMD-II with parameters of IQl 
and those of proximity potential, respectively. 

Fig. 3 The static Coulomb barriers of ^^ii+isij^ 54Q^^208pb^ 32g^230rjj^ i2Q^250p^_ 
The solid and dashed curves denote the results of ImQMD-II with parameters of IQl 
and those of proximity potential, respectively. 

Fig. 4 The probability of fusion reaction gfus{E,b) as a function of impact parameter for 
^°Ca+^sCa at E^,rn. = QOMeV. 

Fig. 5 The probability of fusion reaction gfus{E,b) for ^^Ca+'^'^Ca at -Ec.m. = 
60, 58, 54, 52Me1/. 
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..,.6 The ..io„^xcitat.n fun.io,. ot a^Ca^^T, ,28), b,«T..«T.^ c) a„d d) 



Qg^^90,96 2rj25|, e) and f) 32=34g+89Y[3o|^ and g) and h)28Si,35Cl+922i.j3i| at ener- 



gies near the Coulomb barrier. The open circles denote the experimental data, and 
the filled triangles denote the results of ImQMD-II with parameters of IQl. The 
crosses denote the results of one-dimension WKB approximation 

Fig. 7 The fusion excitation function for '^'^Ca+^^Ca. The stars denote the experimental 
data. The solid circles and triangles denote the results without and with the surface- 
symmetry term taken into account, respectively. 



Fig. 8 The fusion excitation function of ^^^Sn+^^Ni 33|. The solid circles denote the exper- 
imental data. The filled triangles denote the results of ImQMD-II with parameters of 
IQl. The dashed and solid curves denote the results of coupled-channel calculations 
including inelastic excitation (IE) of the projectile and target and IE plus neutron 
transfer (n&IE), respectively. 

Fig. 9 The charge distribution of fragments for (a) Ca+Ca a.t E = 35MeV/nucleon and 
(b) Xe+Sn E = 50 MeV/nud eon. The open circles denote the experimental data, the 
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stars and solid circles denote the results of ImQMD-II model without and with the Qt 
term taken into account, respectively. The calculation results are averaged over the 
impact parameter range of 6 = 1 — 3/m. 

Table . 1 parameters used in the QMD model and the corresponding values obtained from 
various Skyrme parametrizations. 

Table . 2 The parameters of IQl. 

Table . 3 The binding energies and root-mean-square charge radii of a series ground state 
nuclei calculated by the ImQMD-II model with IQl interaction. 
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